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ABSTRACT 

We  present  an  0(n  log2n)  algorithm  for  finding  the  shortest  vertical 
segment  ("watchtower")  uv  that  has  to  be  erected  on  a  polyhedral  terrain  S 
with  n  faces,  so  that  its  top  endpoint  can  see  the  entire  terrain  S.  An  exten- 
sion of  our  technique  yields  an  0(n  log2n  +  Jt)-time  algorithm  for  calculating 
the  lower  envelope  of  two  polyhedral  terrains,  one  of  which  is  convex,  where 
k  is  the  combinatoial  complexity  of  the  output  envelope. 

The  Watchtower  Problem 

The  problem  studied  in  this  paper  can  be  stated  as  follows.  Given  a  polyhedral  terrain  S 
in  3-space  (i.e.  a  polyhedral  surface  intersecting  each  vertical  line  in  exactly  one  point),  find 
the  smallest  length  A  of  a  vertical  segment  uv  (a  "watchtower")  whose  lower  endpoint  u  lies 
on  S  and  whose  upper  endpoint  v  can  see  the  entire -terrain  S  (i.e.,  the  segment  connecting  v 
to  any  point  in  S  lies  entirely  above  S).  See  Figure  1  for  an  illustration.  This  is  one  problem 
in  a  collection  of  related  visibility  problems  for  polyhedral  terrains,  which  have  recently  been 
considered  in  [DFPAN],  [CS].  These  problems  have  obvious  applications  in  processing  real- 
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life  topographic  data,  approximated  by  a  polyhedral  terrain,  for  navigation,  surveillance, 
graphic  rendering,  etc. 

We  can  reformulate  the  problem  as  follows.    Let  flt  .  .  .  ,/„  be  the  planar  faces  of  S, 

and  let  iri ir„  denote  the  planes  containing  these  faces.    It  is  easy  to  see  (as  was 

already  argued  in  Cole  and  Sharir  [CS])  that  a  point  v  can  see  the  entire  terrain  S  if  and  only 
if  it  lies  above  each  of  the  planes  ir,.  Thus  the  region  of  points  that  see  the  entire  S  is  simply 

the   intersection   of   the    upper   half-spaces    it i" -n„     bounded    below    by    irj,  .  .  .  ,irn 

respectively.  Denote  this  region  by  K.  K  can  be  easily  calculated  in  time  0(n  log  n)  (see 
Preparata  and  Muller  [PM]).  It  is  a  convex  unbounded  polyhedral  region,  whose  boundary 
L  =  dK  is  a  convex  polyhedral  terrain  having  at  most  n  faces. 

The  problem  can  now  be  formulated  as:  Given  a  polyhedral  terrain  S  with  n  faces,  and 
another  convex  polyhedral  terrain  L  with  n  faces  lying  above  S,  find  the  shortest  vertical  dis- 
tance between  S  and  L,  and  find  points  u€S,  v€L  which  realize  this  distance. 

(Note  that  we  assume  now  that  L  and  S  do  not  intersect.  Cole  and  Sharir  show  how  to 
test  whether  L  and  S  do  intersect,  in  0{n  log  n)  time,  and  also  how  to  find  a  point  in  the 
intersection,  if  it  exists.  In  this  case  the  vertical  distance  between  L  and  S  is  zero,  and  the 
point  in  their  intersection  as  produced  in  [CS]  can  be  taken  as  the  (zero-length)  vertical  seg- 
ment realizing  this  zero  distance.) 

We  present  below  an  algorithm  that  solves  this  problem  in  time  0(n  log2/!).  We  then 
extend  our  technique  to  obtain  an  algorithm  for  the  calculation  of  the  lower  envelope  (i.e. 
pointwise  minimum)  of  two  polehdral  terrains,  one  of  which  is  convex,  in  time 
0(n  log2n  +  k),  where  k  is  the  combinatorial  complexity  of  the  resulting  envelope. 

The  Algorithm 

The  overall  structure  of  our  algorithm  is  quite  simple.  It  is  easily  checked  that  the  verti- 
cal distance  between  L  and  S  is  realized  by  a  segment  uv,  with  «€5,  v€L,  such  that  either 

(a)     u  is  a  vertex  of  J;  or 


(b)  v  is  a  vertex  of  L;  or 

(c)  u  lies  on  an  edge  of  5  and  v  lies  on  an  edge  of  L. 

It  is  easy  to  handle  the  first  two  cases.  We  project  5  and  L  onto  the  xy-plane,  to  obtain 
two  planar  maps  S* ,  L*  respectively,  and  we  preprocess  each  of  them  for  fast  point  location 
(as  in  Kirkpatrick  [Ki],  Edelsbrunner,  Guibas  and  Stolfi  [EGS],  or  Sarnak  and  Tarjan  [ST]). 
Then,  for  each  vertex  u  of  S  we  locate  its  xy-projection  u*  in  L* .  If /is  the  face  of  L  whose 
projection  contains  «*,  then  the  vertical  distance  between  u  and  L  is  the  vertical  distance 
from  u  to  /,  which  can  be  computed  in  constant  additional  time.  Thus,  in  time  0{n  log  n)  we 
can  calculate  the  shortest  vertical  distance  between  a  vertex  of  S  and  L,  and  similarly  we  can 
calculate  the  shortest  vertical  distance  between  a  vertex  of  L  and  S  within  the  same  time 
bound. 

The  difficulty  lies  in  case  (c).  Of  course,  we  could  test  each  pair  of  edges,  one  from  S 
and  one  from  L,  whether  their  xy-projections  intersect,  and  if  so  calculate  their  vertical  dis- 
tance, and  thus  obtain  the  smallest  such  distance  in  0(n2)  time.  Below  we  show  how  to 
improve  this  crude  approach  to  an  efficient  0(n  log2n)  procedure. 

Let  e  =  ab  be  an  edge  of  S.  For  each  Ostrsl  define  u(t)=  (l-/)a  +  tb  €  e,  and  v(r)  = 
the  point  on  L  lying  directly  above  w(f).  Let  Ft(t)=  |«(r)v(r)|.  It  is  easily  checked  that  Fe  is 
a  piecewise  linear  convex  function  of  t  (where  breakpoints  on  Fe  occur  when  v(f)  lies  on  an 
edge  of  L).  It  therefore  has  a  global  minimum  which  is  attained  either  at  a  breakpoint  or  in 
some  subinterval  of  e  between  adjacent  breakpoints.  Moreover,  given  any  point  u(t)  on  e,  if 
we  know  the  corresponding  point  v(»)  and  the  face  of  L  on  which  it  lies,  then  we  can  deter- 
mine, in  constant  additional  time,  on  which  side  of  u  (t)  the  minimum  of  Fe  is  attained,  or 
that  this  minimum  is  attained  at  t  itself.  (Note  also  that  finding  v(f),  given  «(f),  can  be  done 
in  0(log  n)  time,  by  locating  the  projection  u*(t)  of  u(t)  in  the  map  L* .) 

To  elaborate  this  observation,  let  0<f!<r2<  '  ■  •  fm<l  be  the  "breakpoints"  of  Fe,  so 
that  Fe  is  linear  over  each  of  the  subintervals  [0,r,],  [tut2],  .  .  .  ,[rm,l].  Note  that  the  f/'s 
correspond  exactly  to  points  v(^)  that  lie  on  edges  of  L.    Thus,  by  a  binary  search  over  the 
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sequence 

t\,  .  .  .  ,tm  we  can  determine  the  minimum  of  Fe,  which  must  be  attained  at  one  of  these 
points.  Since  m=0(n),  this  procedure  will  take  overall  0{\o£n)  time.  The  problem,  how- 
ever,  is  that  we  do  not  know  in  advance  the  breakpoints  tit  nor  can  we  afford  to  calculate  all 
of  them,  because  that  would  lead  to  an  overall  quadratic  complexity. 

To  overcome  this  difficulty,  we  proceed  as  follows.  Let  c\,  .  .  .  ,cp  be  the  vertices  of  L, 
and  let  x\,  .  .  .  ,xp  be  their  x-coordinates,  sorted  in  increasing  order.  Let  e  =  ab  be  an  edge 
of  S,  and  let  a^<b\  be  the  x-coordinates  of  its  endpoints.  In  O(logn)  time  we  can  find 
indices  i,j  such  that 

x,_i<a,=s.r,=5xJs*1<x/  +  1  . 
We  then  run  a  binary  search  on  the  points  x,,  .  .  .  ,Xj.  At  each  tested  point  x*  we  find  the 

point  u  (f)  on  e  whose  x-coordinate  is  x*,  and  then  apply  the  above  procedure  to  determine,  in 

O(log  n)  time,  whether  the  minimum  of  Fe  is  attained  to  the  left  or  to  the  right  of  (or  at) 

u(t).  Thus,  in  O(log  n)  such  steps  (of  overall  cost  0(log2n))  we  can  obtain  two  adjacent 

values  Xk<Xk+i  such  that  the  minimum  of  Fe  is  attained  at  a  point  u(r)  whose  x-coordinate 

lies  between  x*  and  xl  +  1. 

Consider  the  projection  L*  of  L  on  the  xy-plane.  The  slab  a  :xt<x<xt  +  i  contains  no 
vertex  of  L  * ,  and  so  its  intersection  with  L  *  consists  of  a  sequence  of  (portions  of)  projected 
edges  g* ,  .  .  .  ,g*  oi  L* ,  each  extending  from  the  left  boundary  to  the  right  boundary  of  the 
slab  or.    Suppose  g*  ,  .  .  .  ,g*  are  sorted  in  increasing  y  order;  see  Figure  2. 

Let  w  and  z  be  the  intersections  of  the  projection  e*  of  our  segment  e  with  the  slab 
boundaries  x=xt,  x  =  xt+i  respectively.  Suppose  w  lies  between  gf-i  and  g*  and  z  lies 
between  g*  and  gj+l  (where  we  assume  without  loss  of  generality  that  z'<v').  Then  it  is  easy 
to  check  that  within  a,  e*  intersects  exactly  the  edges  g* ,g*+\,  .  .  .  ,gf .  Hence  again  we  can 
apply  the  binary  search  technique  described  above,  to  find  the  intersection  at  which  the 
minimum  of  Fe  is  attained. 
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To  facilitate  the  availability  of  this  information,  we  use  the  point  location  technique  of 
Sarnak  and  Tarjan  [ST].  In  their  approach,  the  edges  of  L*  crossing  a  slab  cr  as  above  are 
stored  in  some  binary  search  tree  (in  which  they  are  ordered  by  increasing  y-order).  Using 
persistent  search  trees,  they  are  able  to  store  all  these  search  trees,  one  per  slab,  in  overall 
O(n)  space,  without  degrading  the  logarithmic  search  time  within  each  tree. 

For  us  this  is  a  useful  structure  that  we  can  exploit  as  follows.  We  search  w  and  z  in  the 
tree  Ta  corresponding  to  the  slab  a,  thereby  obtaining  (in  0(log  n)  time)  the  indices  /  and  j 
defined  above.  The  segments  g* ,  .  .  .  ,g*  are  leaves  of  O(log  n)  disjoint  subtrees  of  T„, 
whose  roots  are  children  of  the  nodes  along  the  search  paths  to  w  and  to  z;  see  Figure  3. 

We  can  therefore  check  each  of  these  0(log  n)  trees,  to  find  which  one  contains  the 
segment  g£  whose  intersection  with  e*  gives  us  the  point  at  which  Fe  is  minimized.  (For 
example,  we  can  test  e*  against  the  first  and  last  segments  in  each  subtree.)  Having  found 
this  subtree,  we  can  proceed  with  binary  search  through  it  to  find  the  desired  segment  g% . 
All  this  can  be  done  in  0(log2n)  time  -  O(logn)  binary  search  steps,  each  requiring 
0(log  n)  time. 

Our  overall  algorithm  simply  applies  this  search  to  each  edge  e  of  S,  and  chooses  the 
one  yielding  the  smallest  minimum  vertical  distance.  We  thus  obtain 

Theorem  1:  Given  a  polyhedral  terrain  S  having  n  faces,  the  shortest  vertical  watchtower  on 
S  whose  top  endpoint  can  see  the  entire  terrain  5,  can  be  calculated  in  time  0(n  log2n),  and 
in  linear  space. 

The  Lower  Envelope  of  Two  Polyhedral  Terrains 

Our  method  can  be  generalized  to  solve  the  following  problem.  Given  two  polyhedral 
terrains  L  and  S  having  a  total  of  n  faces,  so  that  L  is  convex,  calculate  the  lower  envelope  M 
of  L  and  S  (i.e.,  viewing  L  and  5  as  the  graphs  of  two  corresponding  functions  z=fL{x,y), 
i=fs(x,y),  M  is  the  graph  of  the  pointwise  minimum  of  fL  and  fs).  Clearly,  M  is  also  a 
polyhedral  terrain,  but  its  combinatorial  complexity  could  be  as  high  as  H(n2),  as  in  the  case 
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shown  in  Figure  4.  Our  goal  is  to  calculate  M  by  an  output-sensitive  algorithm,  whose  time 
complexity  depends  on  the  combinatorial  size  of  M,  and  we  achieve  this  goal  by  using  a 
modified  version  of  our  watchtower  algorithm.  The  lower  envelope  problem  studied  here  is  a 
special  case  of  the  more  general  problem,  in  which  L  and  S  are  arbitrary  polyhedral  terrains; 
we  do  not  know  of  a  similar  output-sensitive  algorithm  for  the  genral  case. 

Consider  first  the  problem  of  calculating  the  vertices  of  M .  Clearly,  each  vertex  is  either 
(i)      a  vertex  of  S  lying  below  I; 
(ii)     a  vertex  of  L  lying  below  S; 

(iii)    an  intersection  of  an  edge  of  S  with  a  face  of  L;  or 
(iv)    an  intersection  of  an  edge  of  L  with  a  face  of  S. 

The  first  two  types  of  vertices  can  be  found  in  time  0(n  log  n),  using  planar  point  loca- 
tion as  above.  The  third  type  of  vertices  is  found  by  noting  that  each  edge  of  S  cuts  L  in  at 
most  two  points,  and  that  these  points  can  be  found  in  0(log2n)  time,  using  the  same  binary 
search  technique  described  above. 

To  find  vertices  of  the  fourth  type  is  a  bit  trickier,  because  an  edge  of  L  might  intersect 
many  faces  of  S  (e.g.  as  shown  in  Figure  4).  We  overcome  this  problem  by  first  collecting  all 
faces  of  S  which  intersect  L.  Without  loss  of  generality,  assume  that  both  I  and  S  are  tri- 
angulated. Then  a  face  A  of  S  intersects  L  if  and  only  if  either  an  edge  of  A  intersects  L,  or 
all  edges  of  A  lie  below  L  and  a  vertex  of  L  lies  below  A.  This  condition  can  be  easily 
checked  as  a  by-product  of  our  search  for  vertices  of  type  (i)-(iii)  at  no  additional  overhead. 

Next,  for  each  face  A  of  5  that  intersects  L,  we  apply  a  "local  search"  to  find  this  inter- 
section. Let  us  denote  by  A/,  the  (convex)  intersection  of  A  with  the  convex  half-space  K 
lying  above  L,  and  let  A£  be  the  projection  of  this  intersection  on  the  j^-plane.  Clearly  the 
intersection  ADL  projects  to  the  portion  t  of  the  boundary  of  A£  that  lies  in  the  interior  of 
the  projection  A*  of  A.  Note  also  that  if  no  edge  of  A  intersects  L  then  t  is  a  closed  convex 
curve  lying  entirely  within  the  interior  of  A*;  otherwise  t  will  consist  of  several  connected 
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arcs,  each  starting  and  ending  on  the  boundary  of  A*. 

We  obtain  t  by  finding  a  starting  point  along  each  of  its  connected  portions  and  then 
trace  t  from  these  points.  If  t  meets  the  boundary  of  A*,  we  find  these  starting  points  by 
intersecting  the  edges  of  A  with  L  (in  0(log2n)  time  per  edge).  Otherwise,  we  take  a  vertex  v 
of  L  lying  below  A,  find  (in  0(log  n)  time)  a  point  u  on  A  lying  directly  above  v,  take  any 
ray  p  emerging  from  u  along  A,  and  find  its  intersection  with  L,  whose  projection  on  the  ry- 
plane  will  necessarily  be  a  point  on  t,  in  time  0(log2n). 

Let  ?*  be  such  a  starting  point  on  t.  We  trace  t  from  z*  in  counterclockwise  direction, 
as  follows.  We  find  the  point  z  on  A  D  L  which  projects  to  z* ,  and  the  face  f  of  L  containing 
z.  We  then  follow  /  along  the  straight  edge  /  ("I  A  in  counterclockwise  direction,  until  its 
intersection  with  the  boundary  of  /  (note  that  this  intersection  projects  to  a  vertex  of  t).  We 
then  repeat  this  edge  tracing  along  the  new  face  of  L  adjacent  to  /  at  this  intersection,  and 
continue  until  either  we  cross  an  edge  of  A  or  we  come  back  to  the  starting  point  z.  Clearly 
this  tracing  can  be  accomplished  in  time  proportional  to  the  number  of  edges  and  vertices  of 

T. 

We  omit  further  details  concerning  the  completion  of  the  calculation  of  M,  and  leave  it 
to  the  reader  to  verify  the  summarizing  theorem 

Theorem  2:  The  lower  envelope  of  two  polyhedral  terrains  L  and  S,  having  a  total  of  n  (tri- 
angular) faces,  so  that  L  is  convex,  can  be  computed  in  time  0(n  log2n  +  k),  where  k  is  the 
combinatorial  complexity  of  the  envelope,  and  in  0(n  +k)  space. 

Discussion  and  Open  Problems 

In  this  section  we  conclude  with  some  comments  and  some  open  problems. 

(1)  The  watchtower  problem  solved  here  is  a  relatively  easy  special  case  of  more  general 
problems  in  which  we  wish  to  find  a  collection  of  short  watchtowers  whose  tops  can  collec- 
tively see  the  entire  terrain.  Many  variants  of  this  problem  suggest  themselves,  with  varying 
criteria  for  optimality  (e.g.  minimum  number  of  towers,  smallest  total  length  of  all  towers, 
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etc.)-  One  suspects  that  many  of  these  problems  are  NP-hard,  as  was  shown  in  [CS]  for  a 
simple  variant  of  this  general  problem,  in  which  we  wish  to  find  the  smallest  number  of 
points  lying  on  a  terrain  S  which  collectively  can  see  the  entire  S. 

(2)  Dobkin  and  Kirkpatrick  [DK]  have  developed  several  algorithms  to  detect  intersection  or 
proximity  between  convex  polyhedra  in  three  dimensions.  However,  their  data  structures  do 
not  seem  to  be  appropriate  for  our  problem,  either  because  the  query  time  needed  to  access 
these  structures  is  too  high,  or  because  they  do  not  address  the  problem  of  finding  the  shor- 
test vertical  distance  between  two  such  polyhedra,  which  is  what  we  require. 

(3)  Re-examining  the  abstract  problem  that  has  actually  been  solved  in  this  paper,  namely  the 
problem  of  finding  the  shortest  vertical  distance  between  two  polyhedral  terrains  when  one  of 
them  is  convex,  an  obvious  generalization  would  be  to  ask  the  same  question  for  an  arbitrary 
pair  of  polyhedral  terrains.  The  generalized  problem  is  considerably  more  difficult  than  the 
one  studied  in  this  paper.  It  follows  from  the  results  in  [ChS]  that  it  can  be  solved  in  subqua- 
dratic  time  (albeit  very  close  to  quadratic),  and  it  remains  an  open  problem  to  find  a  substan- 
tially better  solution  (the  preliminary  results  in  [CEGPS]  might  be  useful  for  this  problem). 

(4)  A  related  open  problem  is  to  obtain  an  output-sensitive  algorithm  for  calculating  the 
lower  envelope  of  two  arbitrary  polyhedral  terrains. 

(5)  Finally,  the  abstract  problem  studied  here  might  be  too  general  for  the  specific 
watchtower  problem  that  has  originated  it,  because  we  are  not  taking  advantage  of  the  special 
relationship  between  the  two  surfaces  S  and  L.  Can  this  relationship  be  exploited  to  obtain  a 
faster  algorithm? 
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Figure  Is  The  shortest  watchtower  problem 
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Figure  2:  The  structure  of  a  slab  <r 


-  12 


Figure  3:  Searching  through  the  tree  Tc 
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Figure  4:  The  lower  envelope  of  one  convex  and  one  arbitrary  polyhedral  terrains 
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